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Dynamics of neutral atoms in nonuniform magnetic fields, typical of quadrupole magnetic traps, 
is considered by applying an accurate method for solving nonlinear systems of differential equations. 
This method is more general than the adiabatic approximation and, thus, permits to check the 
limits of the latter and also to analyze nonadiabatic regimes of motion. An unusual nonadiabatic 
regime is found when atoms are confined from one side of the z -axis but are not confined from 
another side. The lifetime of atoms in a trap in this semi-confining regime can be sufficiently long for 
accomplishing experiments with a cloud of such atoms. At low temperature, the cloud is ellipsoidal 
being stretched in the axial direction and moving along the z -axis. The possibility of employing 
the semi-confining regime for studying the relative motion of one component through another, in a 
binary mixture of gases is discussed. 

03.75.-b, 02.30.Hq 



I. INTRODUCTION 

The motion of neutral atoms in nonuniform magnetic fields is important to study for several applications, in 
particular, for better understanding the mechanics of confinement in quadrupole magnetic traps, such as the Ioffc- 
Pritchard traps with a static bottle field [1-3] or dynamic traps with a rotating bias field [4,5]. This has become 
especially interesting after the experimental observation of Bose condensation in very cold gases of Rubidium [6,7], 
Lithium [8], and Sodium [9-11]. There exists extensive literature considering statistical properties of confined Bose 
systems, using various approaches, e.g. the quasiclassical density-of-state approximation [12,13], the Gross-Pitaevskii 
equation [14-16], the Monte Carlo density-matrix calculations [17], the Thomas-Fermi approximation [18,19], the 
Bogolubov approximation [20], and the gas approximation in the frame of the Gibbs ensembles [21,22]. Statistical 
properties of the weakly interacting Fermi gas confined in a potential well have also been studied in the Thomas-Fermi 
approximation [23,24]. 

The aim of this paper is to consider not statistics but dynamics of atoms in nonuniform magnetic fields. When one 
is interested in the behaviour of confined atoms, one deals with their stationary motion. Stationary regimes can be 
described by the adiabatic approximation. When the confining potential is harmonic, then the dynamics of atoms is 
given by simple harmonic oscillations. In general, the confining potential is not necessarily harmonic. For instance, 
the first demonstrated magnetic trap [25] used a quadrupole field (with zero magnetic field at the center) which gave 
rise to a linear potential. In any case, atomic motion in a strictly confining potential can be described by the adiabatic 
approximation. Such adiabatic motion in various magnetic traps has been analysed in Refs. [26-29]. 

A more general consideration of atomic motion, without using the adiabatic approximation, is meaningful for several 
reasons: First of all, a more general approach makes it possible to understand the limits of the adiabatic approximation. 
Second, studying other, nonadiabatic, regimes of motion permits to explain more profoundly the physics of atoms 
inside magnetic traps, as far as in these not all atoms are confined. Knowing better different, including nonadiabatic, 
regimes of atomic motion may, possibly, give a hint on how to improve confining characteristics of magnetic traps. 

One more reason is related to the recently reported experiments on the simultaneous trapping of two different atomic 
species, two isotopes of rubidium [30] and sodium and potassium [31]. These experiments are a starting point for a 
new series of studies of ultracold matter. The variety of effects that can be observed in mixtures are incomparably 
richer than in one-component gases. This concerns even equilibrium mixtures [32]. Much more interesting features 
appear when one of the components can move through another. For instance, in a binary mixture with such a relative 
motion the effect of conical stratification [33] can occur. This effect happens when one of the components moves in one 
direction through another component. Then the instability can develop inside a cone with the axis along the relative 
velocity; as a result of this instability, the components stratify in space. Note that for this effect the one-directional 
relative motion is necessary, but not relative oscillations or collisions of sloshing clouds. In the presence of a relative 
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macroscopic motion some unusual manifestations of the Dopplcr effect [34] may also arise. But can such a strange 
regime exist when the atoms of one kind move in one direction being confined from another? In addition, it is desirable 
that this one-directional escape from a trap would not be too fast in order to be able to accomplish measurements. 
These requirements look too severe to allow the existence of such a semi-confining regime. However, in what follows it 
will be shown that this semi-confining regime does exist. Certainly, it is nonadiabatic and, even more, nonpotential, 
thus cannot be described in the frame of the adiabatic approximation. 



II. EVOLUTION EQUATIONS 

Since the aim of this paper is to present an accurate solution of evolution equations for atoms in nonuniform 
magnetic fields, it is reasonable, first of all, to pay attention to the accurate formulation of the equations themselves. 
The quantum Hamiltonian of a system of TV neutral atoms, each with a mass m and magnetic moment /i , is 

N I ^ 2 \ N 

»=1 \ / i& 

where Pi is a momentum operator; Si , a spin operator; Bi=B ( r i,t) is the total magnetic field acting on an 
i -atom; and $y = <f>(| r j — r 3 ■ |) is an interaction potential. The magnetic moment ^ is the product of the 

Bohr magneton and the hyperfine g -factor [35]. The wave function of the system, $ = \^ s { r i, r 2, ■ ■ ■ , r at, t)] is 
a column in spin space. The quantum-mechanical average of an operator from the algebra of observables, A , is the 
scalar product 

(A) = (*,AV). (2) 

Using the Schrodinger equation ihd^f jdt = , with the Hamiltonian (1), it is straightforward to get the evolution 
equations for the average position of an atom, 



its average momentum 



d 
dt 



(P l ) = M(V, {SiBi)) - < v * ( 4 ) 



and the average spin 



^(S^^xBi). (5) 



These equations may be simplified with a mean-field approximation 

(S?Bf) = (S? )(Bf), (6) 

where a,(i — x,y,z , valid for fields slowly varying in space [35]. Condition (6) is also called the semi-classical 
approximation and is usually supplemented by another approximate equation 

(B(ri,t))=B((u),t) (7) 

which again assumes a slow variance of magnetic field in space. When the magnetic field is a linear function of real- 
space coordinates, then Eq.(7) is not an approximation but an exact relation. This concerns quadrupolc magnetic 
fields which in what follows we shall deal with. Hence, the sole approximation we need is the mean-field one Eq.(6). 

Under conditions (6) and (7), equations (3)-(5) acquire the same form for all indices i , which permits us to simplify 
the notation by introducing 

r= (r.) = {x,y,z}, 
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and the average interatomic force 



v=^{Pi) = {v x ,Vy,v z }, (8) 



S (Si) ~ {S X , Sy, S Z }, 



N 

y= - e & (9) 

Then Eqs.(3)-(5) can be reduced to the system of equations 

(I r -> 

nr =v > 

d 4 = ^CsB) + l (10) 

at m m 

which will be the main object of our consideration. 

The derivation of (10), though simple enough, contains an important point which is worth emphasizing. The basic 
approximation (6) supposes that the system fields slowly vary in real space, so that the spin and real-space degrees of 
freedom can be separated, which in the quantum-mechanical language means that the wave function can be factorized 
into a product of spin and real-space wave functions [35] . This quantum-mechanical separation of variables, as will 
be shown in what follows, is closely related to their dynamical separation. 

The evolution equations (10) are to be supplemented by the initial conditions 

r (0) = ?o= {xo,Vo, z }, 
v (0)=v ={v x o ,vlv* Q }, (11) 

5(0) =So={SS,Sg,S5}. 

In specifying the form of the magnetic field, let us take it as in the experiments [5-7] with dynamical quadrupolc 
traps. Then the total magnetic field 



is the sum of the quadrupole field 



and the rotating bias field 



B=Bi [r)+B2 (t) (12) 



Bi (r)=B[(r -3z e x ) (13) 



B2 (t) = B 2 (e x cosojt+ e y sinuit), (14) 



where e a is a unit vector for a = x, y, z . 
The characteristic length 



L-% 05) 



of the quadrupole-field nonuniformity corresponds to the radius of the field zero in the radial direction. This length 
defines approximately the upper limit for the radius of a trapped atomic cloud. Keeping this in mind, it is convenient 
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to pass to dimensionless space variables measuring the components of the Cartesian vector r = {x, y, z} in units of 
L . Then we can make profit from the inequality 

| 7 | < 1. (16) 

To return to dimensional space variables, we have to put r^rL. 
We introduce the characteristic frequencies 

^ = ' W2 = — (17) 

of atomic and spin motions, respectively, and the collision rate 7 defined by the ratio 

in which £ is treated as a stochastic variable representing the interactions of atoms through their random collisions. 

A detailed definition of the variable £ will be given in Sec.V and in Appendix. Then from Eq.(lO) we obtain the 
evolution equations for the space variable 



d 2 r 
dt 2 



and for the spin variable 



= cuj (S x 7 X +S y e„ -2S Z e z ) + 7 ? (19) 



W = ^ A ^' (20) 

where the matrix A = [A a p] , with a, /3 = 1,2,3, consists of the elements 

An = A 22 = A 33 = 0, 

A u = -A21 = -2z, 

A13 = -A31 = -y- sinivt, 

A23, = —A32 = x + cos ivt. 

We are concerned about an accurate solution of the system (19) and (20) of nonlinear differential equations, without 
using the adiabatic approximation. At this point, to distinguish what is what, it is useful to say several words about 
the adiabatic approximation. Fortunately, this will not take too much space, since the latter approximation is rather 
trivial. First, one assumes that the dynamical process is close to its stationary state, so that it is admissible to put 

d S /dt = . Then, from the third equation in (10) it follows that S x B= . This means that the spin is 
aligned along B , which can be written as S • B— S\ B | • Thus, one excludes the spin motion saying that spin 
adiabatically follows the magnetic field. Using this, for the atomic space variable r one gets the Newton equation 
with the adiabatic force 

F a =S(j,V I B |. 

Assume that the field rotates much faster than the mechanical oscillations of atoms, but not so fast as to induce 
transitions in the Zeeman substructure. This means that u>i <C u> <C <jj 2 ■ Then the adiabatic force can be averaged 
over the period 2n/oj of the rotating bias field. Averaging this force and using Eq.(16), one gets 

(F a )t ^ -^Swf(x ~e x +y ~e y +8z ~e z ), 
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which immediately yields the adiabatic potential 

U a ~ ~Su 2 (x 2 + y 2 + 8z 2 ) + const . 

This is a harmonic, although anisotropic, potential. The motion of atoms in such a potential is given by simple 
harmonic oscillations, if S < . When S = , atoms fly away ballistically and if S > , they escape by the 
exponential law. In both latter cases atoms escape from the trap in all directions. The ballistic flying away is isotropic. 
The exponential escape, because of the anisotropy of the adiabatic potential, is anisotropic: atoms escape faster along 
the axial direction than in the radial one; but anyway the symmetry with respect to the inversion r^> — r is 
preserved. 

III. SCALE SEPARATION 

Return to the general equations (19) and (20). Written in the standard form, they compose a nonlinear dynamical 
system of the ninth order, that is, a system of nine nonlinear differential equations. It seems that it is impossible to solve 
this complicated system of equations without invoking a rough approximation like the adiabatic one. Nevertheless, 
these equations can be solved using the method of scale separation [36,37]. The mathematical foundation of this 
approach is based on the Krylov-Bogolubov averaging method [38] and the Poincare theory of generalized asymptotic 
expansions [39]. The method of scale separation was successfully applied to several intricate problems, such as the 
origin of self-organized spin superradiance in nuclear magnets [36], coherent radiation regimes of spin masers [40,41], 
and fast polarization reversal in proton targets used for studying beam scattering [42] . The accuracy of this approach 
has been confirmed by good agreement of its solutions with experimental data [43,44] and with computer simulations 
[45,46]. 

The first step of the method of scale separation [36,37] is to classify the functional variables of the problem, separating 
time scales. Fortunately, in the majority of interesting physical problems it is possible to separate relatively slow from 
relatively fast variables. In our case, such a separation is naturally related to the mean-field condition (6) used 
for deriving the evolution equations (19) and (20). As has been already discussed above, condition (6) assumes 
that the real-space nonuniformity in the system is, in some sense, small. Now we ascribe an exact meaning to 
this phrase concretizing in what sense the nonuniformity is small. The nonuniformity in the system is connected 
with the quadrupole field (13) and pair interactions in (9), while the rotating bias field (14) is spatially uniform. 
Therefore, what we need is to compare the characteristic parameters related to the corresponding fields. Among these 
characteristic parameters we have u\ and u 2 in Eq.(17), and 7 in Eq.(18). The nonuniformity is weak if the 
characteristic parameters corresponding to nonuniform fields are small as compared to that of a uniform field. The 
latter means nothing but the validity of the inequalities 

lui <C U2, 7 < W2. (21) 

The meaning of the first inequality in (21) is quite evident, indicating that the frequency of mechanical oscillations 
of atoms is much smaller than that of spin fluctuations. The second inequality is also very natural, as far as the 
collision rate is usually much smaller than L02 . If 7 would be comparable with u>2 this would imply atomic collision 
are causing radio-frequency transitions between magnetic sublevels. With Eq.(21) in mind, looking at the evolution 
equations (19) and (20), we notice at once that the variable r is to be treated as slow, compared to the fast variable 

S ■ 

At this point, it is worth emphasizing how naturally the separation of functional variables into slow and fast, 
with respect to time, is connected with the character of nonuniformity in real space. This is why spending some 
time for remembering the derivation of Eqs.(10) was not in vain, but, on the contrary, is important for stressing the 
self-consistency of the approximations used. 

Following further the method of scale separation [36,37], we need to solve the equations for fast variables, with slow 
variables being kept as quasi-integrals of motion. The evolution equation for fast variables is Eq.(20) for spin. This 
equation for an arbitrary given antisymmetric matrix A , with elements = — Aji , can be solved exactly. This 

means that we are able to present an exact solution for S for any given external fields. Because of the significance 
of such s solution, we write it down explicitly. 
First, we solve the eigenproblem 

A~bi=ai^ u |T 4 | 2 = 1, (22) 
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in which A is an antisymmetric matrix and i = 1, 2, 3 . The solution is straightforward giving the eigenvalues 

cui = ia, «2 = —ia, a 3 = 0, (23) 

with 



^12 + ^13 + ^23- 



The eigenvectors are 



6 4 = 



1 



with the normalization constant 



(etiA 13 + A 12 A 23 ) e x + (aiA 23 - A 12 A 13 ) e y + (of + A\ 2 ) e 



(24) 



C 1 = (|o i | a -A; 9 ) a + (|a i | a + i4? a ) (Aj 3 + Al 3 ). 



It can be checked straightaway that the vectors from Eq.(24) form an orthonormal basis and satisfy the properties 

(25) 



bibj=5ij, 61=62, 63=63- 
Therefore the general solution of Eq.(20) can be written as a linear combination 
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S (t)=J2 a * Si (t) 



i=i 



of particular solutions 



S t (t)=bi {t)exp{<pi(t)}, 



(26) 



(27) 



in which hi are given by Eq.(24). The coefficients in Eq.(26) are defined by the initial condition for spin from 
Eq.(ll), which yields 



a t =Sob t (0). 
Substituting Eq.(27) into Eq.(20), we obtain the phase 



<Pi(t) 



uJ 2 ai(t)- 6, (*)- 6, (t) 



dt. 



From Eq.(29), invoking Eq.(25), we find that 

ip\ = -(fx, ip* 2 = -tp 2 , tp 3 = 0. 



(28) 



(29) 



(30) 



Let us accent that, of course, not each system of equations like (20) can be solved exactly. The possibility of 
obtaining here the exact solution (26) is due to the antisymmetry of the matrix A . 



IV. ATOMIC VARIABLES 



At the next step of the method of scale separation [36,37] the solution for fast variables is to be substituted into the 
equations for slow variables with time averaging the right-hand side of the latter equations. As follows from Eqs.(27) 
and (29), the fast spin fluctuations are described by the effective time-dependent frequencies ipi(t)/t . Hence, the 
solution (26) does not have a definite period. Consequently, the time averaging is given by the rule 

(F) t = lim - / F(t)dt. (31) 

Averaging Eq.(26) according to the rule (31) and using Eq.(16), we find 
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(S)t = \ [(1 + x)S$ + yS% - 2zS z ] (x e x +y l y -4z e x ). (32) 

The equation for the guiding centers of atomic variables is obtained by time-averaging the right-hand side of 
Eq.(19) with £ treated as a slow variable. Employing Eq.(32), we have 

d 2 7 - -» 

^=F +7 ^, (33) 

with the force 

F= ^ [(1 + x)S* + yS% - 2zS z ] (x t x +y 7 y +8z e x ). (34) 

This force essentially depends on the initial polarization of spins, as it should be in the general case and in contrast 
to the adiabatic force mentioned at the end of Sec. II. Only for one type of initial polarization, when Sq — S and 
Sq = S§ = , the force (34) reduces to the adiabatic one. This type of spin polarization leads to the stationary 
confined motion for which the adiabatic approximation is admissible. 

Let us analyze another situation when spins are initially polarized along the z -axis, so that 

S%=0, S% = 0, S Z = S. (35) 

Then the force (34) becomes 

F= -SJ(z(r +7z ~e z ). (36) 

As is evident, the force (36) is nonadiabatic and, even more, it is nonpotcntial since there exists no potential U such 

that F would be equal to — V U . 
With the force (36), Eq. (33) becomes 

^ + Su>lz{r +7z e x ) = 7 ? . (37) 

The force 7 £ , according to (18), originates from interatomic interactions in Eq.(9). Because of the isotropy of the 
interaction potential <fr(| r j — r j |) , the force corresponding to these interactions may be presented as an isotropic 
vector, that is, we may write 

?=C(e > x+ e y +e x ). (38) 
Expanding Eq.(37) into components, we get for the x -component 

^ + Su%zx - 7^ (39) 

The equation for the y -component is the same as Eq.(39) with the replacement x — > y . Therefore, we shall consider 
in what follows only one of the radial components. For the axial variable, Eq.(37) gives 

+ 8Su; 2 z 2 = 7 £. (40) 

So, we have to solve the system of nonlinear equations (39) and (40) . 
The general solutions to Eqs.(39) and (40) can be written in the form 

x = xi+x 2 , z = zi + z 2 , (41) 

in which x\ and z\ are the solutions to the corresponding homogeneous equations while x 2 and z 2 are the 
solutions to the nonhomogeneous equations. The homogeneous equations are 



d 2 xi 

~~dW 



+ Suj\z\X\ = 0, 
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d 2 Zl 

dt 2 



+ 8Sujlzf = 



(42) 



with the initial conditions 



x 1 (0)=x , x 1 (0)=v%, 
zi(0) = zo, ii(0) = «o, 



(43) 



where the dot means the time derivative. Writing down the nonhomogeneous equations, simplifying them a little, 
taking into account that the collision rate is typically an order of magnitude smaller than uj\ , yields 



d 2 x 2 
dt 2 

d 2 z 2 
dt 2 



Sui\z\x 2 = 7£ - Suj 2 x\z 2 , 
+ \§S<jj\zyz 2 = 7£; 



(44) 



the initial conditions for Eq.(44) being 



x 2 {0) = 0, i 2 (0)=0, 
z 2 (0) = 0, i 2 (0) = 0. 



(45) 



To solve the system of equations (42) and (44), we have, first, to solve the equation for z\ , than to substitute z\ 
into the equation for x\ , and to use the found x\ and z\ in Eq.(44). 



V. SEMI-CONFINED MOTION 

The system of equations in (42) for X\ and Z\ is a system of two nonlinear differential equations of second order. 
However, these equations can be solved exactly. 

Integrating once the second equation in (42), we get 



dt J 3 lK rn lh 
where z m is an integration constant which can be found from the initial conditions in (43) yielding 

Z m ~ Z + mS(j2 - 



(46) 



(47) 



Since the left-hand side in (46) is nonnegative, this implies that > z\ if S > , and zj^ < z\ if S < , that 
is, z m is the maximal value of z\ for S > and z m is the minimal value of z\ for S < , 



max zi(t), S > 0; 



z m — mmzi(t), S < 0. 



We introduce a function 



for which Eq.(46) transforms to 



p = --s<4z u 



— ) =AV 3 -g 2 V-g 3 , 



(48) 



(49) 



(50) 



where 
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.92 = 0, 



93 



256 J3 o3 6 



(51) 



Eq.(50) is the Weierstrass equation with the Weierstrass invariants in (51) and with the discriminant 

A = g\ - 27gj = -27gj. 

The solution of the Weierstrass equation (50) is the Weierstrass function Vit — to) , where to is an integration 
constant. The Weierstrass function is an elliptic function, that is, a doubly-periodic function which is analytic, except 
at poles, and which has no singularities other than poles in the finite part of the complex plane. All properties of the 
Weierstrass functions are perfectly described in Refs. [47,48]. 

In this way, the solution of Eq.(46), and therefore of Eq.(42), reads 



zi(t) 



T(t-t a ). 



The integration constant t is to be found from the initial condition in (43), which gives 

W(t ) = -4z Sw 2 , 

where we took into account that the Weierstrass function is an even function. When t tends to t , then 



zi(t) ~ - 



7^yi + |(i-io) 4 



(52) 



(53) 



(54) 



The general behaviour of zi(t) is as follows. If S > , then, starting from zq , the value of z\(t) increases to z m 
from (47), after which it decreases, diverging to — oo as t — > to ■ If S < , then z\(f) decreases to its minimal 
value z m , after which it turns to the positive direction, increasing to +oo as t — ► to ■ 
Substituting solution (52) into the first of Eqs.(42), we have the equation 



-^ = - 4 v(t-t o)xi 



(55) 



for the radial motion. This is a Lame equation of degree n = 1/2 , which is defined by the relation n(n + 1) = 3/4 . 
The solution to the Lame equation is given by combinations of Lame functions of different kinds [49] . In the present 
case, the solution to Eq.(55) is 



Xl(t) 



ClP 



t-t 



e: 



-1/2 



t-tr 



(56) 



Here 



E 3 (t) = -V(t) 

is a Lame function of degree 3 , of the first kind. The integration constants c\ and are defined by the initial 
conditions in (43). 

The solution (56) diverges together with (52), as t — > to , by the law 



Xl(t) ~ Ci 



t-t 



-1/2 



+ C 2 



t-t 



3/2 



(57) 



Comparing Eqs.(57) and (54), we see that the divergence along the axial direction is faster than in the radial one; the 
aspect ratio being 



z?(t) 



\t-t \ 3/2 . 



(58) 



This ratio tends to zero, as t — > to ■ Therefore, a cloud of atoms acquires an ellipsoidal shape stretched in the axial 
direction. 
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The behaviour of x\{t) and z\{t) shows that the atoms with an initial polarization S > are confined from the 
side z > but are not confined from the side z < . Vice versa, the atoms with an initial polarization S < 
are confined from the side z < but are not confined from that z > . Thus, as ensemble of atoms, with a given 
polarization, loaded into a trap would move predominantly in one direction either to z < or to z > depending 
on whether the initial polarization is S > or S < , respectively. Such a regime is exactly that semi-confining 
regime we have been looking for. 

However, we need yet to find the solutions of nonhomogeneous equations in (44). To this end, we have to concretize 
the variable £ originating from the averaged pair interactions. For a rarefied system these interactions can be treated 
as random pair collisions. The system is said to be rarefied if the average atomic density p and the scattering length 
a satisfy the inequality pa 3 <C 1 . This is just the case of experiments [5-11] with alkali atoms. Therefore, we can 
consider £ as a random variable. More accurately, £(t) can be modelled by a stochastic field. The latter may be 
specified as white Gaussian noise [50] with the stochastic averages 

«£(*)» =0, (mm))=2D6(t-t% (59) 

where D is a diffusion rate. Then the nonhomogeneous equations in (44) become stochastic differential equations. 
An explicit way of treating the random variable £ is presented in Appendix. 

When solving the nonhomogeneous equations, it is useful to invoke again the idea of scale separation. The time 
variation of solutions to stochastic differential equations reflects the properties of the given stochastic fields. Since in 
our case, the stochastic field is modelled by white noise, which is characterized by sharp time jumps [50], then 

the related functional variables x 2 and zi can be treated as fast compared to x\ and z\ satisfying the equations 
not containing such random fields. Hence, for the stochastic differential equations in (44), the slow variables x\ and 
z\ can be kept as quasi-integrals of motion [36,37]. Then using the method of Laplace transforms, we obtain 

x 2 (t) = f G x (t - t) [ 7 £(t) - SujIx iZ2 (t)] dr (60) 
Jo 

and 

Z2(t)= [ G z {t-T)^{T)dT, (61) 

Jo 

where the initial conditions in (45) have been taken into account, and the transfer functions are 
with the effective frequency 

e = \J Sz\uj\. 

Employing condition (59), we may find the moments of solutions (60) and (61). For instance, 

«s 2 (t)»=0, «* 2 (t)»=0. (62) 
Calculating the second moments, we get the mean-square deviations for the radial random variable, 



«^)» = ^fl-^ 



+ 



T 2 Dix? f sm(4et) 1 
+ 3 g 0Q£ 2 2 | 1 ~ cos ( £t ) cos(4rt) H [cos(et) - cos(4et) - l&et sin(ei)] j , (63) 



and for the axial random variable, 



«*.»> - S 



1 



sin(8e£) 



Set 



(64) 



Note that the collision rate 7 enters here as 7 2 , thence the behaviour of the random variables (63) and (64) does 
not depend on whether the interatomic interactions are repulsive or attractive. 
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When the variables x\ and z\ diverge as in Eqs.(56) and (54), then for the random variables in (63) and (64) we 
obtain 



as t —> t . Since this expansion is governed by the same exponentials, it is practically isotropic, with only a slight 
anisotropy due to different preexponential factors. 

Remember that the general solutions to the equations (39) and (40) have the form of the sums in (41) containing 
both the regular terms x\ and z\ and the random terms x 2 and z 2 ■ The relative contribution of these terms is 
regulated by the relation between the parameters 7, D , and lo\ . If -f 2 D <C uif , then the influence of the random 
terms x 2 and z 2 is negligibly small, and the atomic motion is characterized by the regular terms x\ and z\ . In 
this case we have the semi-confining regime. An ensemble of atoms would form an ellipsoidal cloud moving in one of 
the directions along the z -axis. If j 2 D ^> ojf , then the motion of atoms is governed by the random terms x 2 and 
z 2 ■ In such a case, an ensemble of atoms would form an almost isotropic exponentially expanding cloud. 

Assuming, as usual, that the diffusion rate D is proportional to temperature T , we come to the conclusion that 
the realization of either the regime of the fast exponential expansion or the regime of the slow semi-confined motion 
depends on temperature. At high temperatures the former regime will be realized while at low temperatures, the 
latter. A crossover temperature T c related to the equality -f 2 D = lo\ would correspond to the effective boundary 
between these two regimes. 



In order to impart to the whole consideration a completely realistic flavour and to show the reasonability of all 
inequalities assumed for employing the method of scale separation, let us adduce numerical estimates basing on the 
characteristic quantities typical of the experiments [5-7] with 87 Rb in dynamical quadrupole traps. 

The mass of a Rubidium atom is m — 1.45 x 10~ 22 g and the magnetic moment is fx = 0.45 x 10 -20 erg/G . The 
gradient of the quadrupole field is B[ = 120 G/cm , and the amplitude of the rotating bias field is B 2 = 10 G ; 
the rotation frequency of the latter being lu w 5 x 10 4 s" 1 . The nonuniformity length (15) is L <~ 0.1 cm . 
The characteristic frequency of atomic motion and the Larmor frequency of spins from (17) are uj\ <~ 10 2 s _1 and 
uj 2 <~ 5 x 10 7 s _1 , respectively. The collision rate in (18) can be estimated as 7 ~ hpa/m , which, with the average 
density p~3x 10 12 cm~ 3 and the scattering length a ~ 10~ 6 cto , gives 7 ~ 10 . Thus, the following inequalities 
hold true: 



Hence, the classification of atomic variables as slow and of spin variables as fast, based on the inequalities in (21), is 
correct. 

To estimate the lifetime to of the semi-confined motion described in Sec.V, we have to return to Eq.(50). Integrating 
this over time between t — and some t , we have 





(65) 



VI. NUMERICAL ESTIMATES 



7«wi«w«w 2 . 




From here, taking into account that g 2 — and V(t — t ) 



00 , as t — > to j w e get 




(66) 



Comparing Eqs.(51) and (53), we find for the Weierstrass invariant 
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Then Eq.(66) can be rewritten as 



f z ° dz 

to = r t= _ , (67) 

J — i 



where we assume that S > and 

TO 







J_ /£ 
4wi V S" 



The value of z m is given by Eq. (47). The average kinetic energy )^m,{v^) 2 L 2 can be expressed through temperature 
as \ksT . Therefore Eq.(47) acquires the form 

z 3 m = z 3 +(, (68) 

in which 

3T _ mSwlL 2 

C= 167V ° = A* ' 

Finally, for the lifetime (67) we obtain 

-7?TTk- (69) 

For the quantities considered, we have tq w 10~ 2 s and T ~ 10~ 4 if . The parameter £ depends on temperature. 
In the interval of temperatures between 1 nK and 1 mK , it changes from 10~ 6 to 1 . The quantity zq, according 
to Eq.(43), is the initial coordinate zi(0). It is clear that a cloud of trapped atoms should have a distribution of z 
which depends on the temperature and shape of the confining potential. The average value of z n can be interpreted 
as the initial location of the center of an atomic cloud. This value can be different for different experiments. Keeping 
this in mind, we consider below several values of zq. The integral (69) was calculated numerically for different initial 
positions zq . The results are presented in Table 1. As is seen, for the temperatures T w 10~ 5 K the time to 
becomes of order 0.1 s , and for Twl nK it reaches to w 0.3 s . Moreover, diminishing the temperature lower, 
it is possible to make to arbitrary large, since to — > oo when £ — ► . Notice that the lifetime to does not change 
much when varying z m the interval —0.1 < z a < 0.1. 

It is possible to pose the question: how well does the lifetime t characterize the real escape time of an atom from 
a trap? Looking at Eq.(54), we see that z\{t) diverges as t —* to , while, according to condition (16), the actual 
escape of an atom from a trap occurs at the time t\ , when z(t\) ~ 1 . The relation between the times to an d 
t\ , as follows from Eq.(54), is t\ — to ~ w r 1 . This, with the given lo\ ~ 10 2 s _1 , makes t\ — to ~ 10 _2 s . So, 
if to > 0.1 s , then t\ ~ to . Therefore, the time to really plays the role of the average lifetime of atoms in a trap 
during the regime of semi-confined motion. 

If we are not satisfied by the simple estimates for the characteristic time t\ , we can calculate it exactly from the 
equation of motion. The procedure is the same as that for calculating to . The results of this calculation confirm 
that t\ is very close to the time to , being smaller by about 0.02 s . Because of the mutual closeness of these times, 
we do not repeat for t\ the whole table as for t but, for cogency, we present in Table 2 the values of t\ for the 
particles with the initial location at the centre of the trap. These values are found from the formula 

tl = To / -ITT?- (70) 
Jo v z + c 

Similarly to to , the time t\ also becomes arbitrary large, t\ — > oo , as £ — > , i.e. when temperature decreases. 

Treating the interactions of atoms as random pair collisions is admissible if the atomic system is rarefied, so that 
pa 3 <C 1 . For the case considered, p <~ 10 12 cto~ 3 and a ~ 10~ 6 cm , yielding pa 3 ~ 10~ 6 . Thus, the assumed 
inequality is well satisfied. 

As is argued at the end of Sec.V, at high temperatures there exists a regime of fast exponential expansion of an 
almost isotropic cloud. At low temperatures, the regime changes to the semi-confined motion of a slowly moving 
ellipsoidal cloud. The crossover temperature T c can be defined by the equality -f 2 D = uo 3 . The latter, with the 
diffusion rate D ~ ksT/h , gives 

T ~ 

k B l 2 ' 
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Substituting here ui\ <~ 10 2 s _1 and 7 <~ 10 s _1 , we have T c ~ 10 _7 .K" , which is close to the temperature of the 
Bose condensation observed in experiments [6,7]. Consequently, the semi-confined regime can also be realized under 
similar conditions. 

Keeping in mind the future possibility of applying statistical methods to describe an ensemble of atoms in the 
semi-confining regime, we have to understand whether the local equilibrium can be established in this case. The time 
of local equilibrium is related to the collision rate as r/ oc ~ 7 -1 . Hence for 7 ~ 10 s , we have r/ oc ~ 0.1 s . The 
local equilibrium develops if the lifetime of the cloud, to , is longer than the local-equilibrium time. If to ~ 0.3 s , 
then the local equilibrium can be achieved. Of course, the global equilibrium for a semi confined motion cannot exist. 
Recall that the corresponding force (36) is not confining for a semi-axis of z , and furthermore it is not potential. 
But the possibility of local equilibrium implies that a statistical description for such a regime can be done, e.g., with 
the help of hydrodynamic equations. If the local equilibrium is absent, one has to use kinetic equations. 



The method of scale separation is applied to the dynamics of neutral atoms in nonuniform magnetic fields typical of 
quadrupole magnetic traps. We concentrate our attention on the case of a dynamic quadrupole trap with a rotating 
bias field, as in experiments [5-7]. For the initial spin polarization Sq < , the motion of atoms is confined and well 
described by the adiabatic approximation. 

For the atoms whose spins are initially polarized along the z -axis a novel unusual regime appears where the 
motion is confined in a half-space: when Sq > , the motion is confined from the side z > , and if Sq < , 
it is confined from the side z < . This semi-confined motion is nonadiabatic and nonpotential. An ensemble of 
atoms in a semi-confining regime forms an ellipsoidal cloud stretched in the axial direction and slowly moving along 
the z -axis to the nonconfining side. 

The semi-confining regime can be used for studying mixtures with a relative motion of components. For example, 
one component can have initial spin polarization Sq < , thus, being confined, while another component, with initial 
spin polarization along the z -axis, can be semi-confined, moving through the first. Or one can load into a trap 
two components both being polarized along the z -axis, but one with Sq > while another with < . Then 
both such species will be in the semi-confining regime but moving in opposite directions, one to the negative z - and 
another to positive z -direction. The variant with three components, one of which is confined with the two others 
semi-confined moving in opposite directions, could also be realized. Mixtures with the relative motion of components 
can display a number of interesting features, e.g., the effect of conical stratification [33]. 

The semi-confining regime can also be employed for separating the species with different initial polarizations. For 
instance, if one loads into a trap a mixture of two components, one having the spin polarization up and another down, 
then the trap will act as a separator moving the first component down and the second one up, thus separating them. 

One more application of the semi-confining regime could be for realizing atom lasers, for which a directed motion 
of atoms is necessary. 
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Appendix. Random Variables 

The variable £(t) describing random pair collisions of atoms has been treated as a stochastic variable. Consequently, 
Eqs.(44) are stochastic differential equations. To find their solution, it was necessary to define the stochastic averages 
(59). The explicit definition of the random variable £(t) and of the corresponding stochastic averages can be done in 
the following way. 

Assume that there is a set {<p n (t)} of functions ip n (t) ennumerated by a multi-index n. Let this set be complete 
and orthonormal, 



VII. CONCLUSION 




n 
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Then a function can be represented as an expansion 

n 

Each coefficient is considered as a random variable with a probability distribution p(£ n )- For concreteness, we may 
think of p n (£n) as a Gaussian distribution. 

The stochastic averaging for a function F(£(t)) of the stochastic variable is defined as the functional integral 

<™*))»= [ Fm)UMZn)dt n . 

J n 

If the random coefficient £„ is complex, then d£ n = d(Re£ n )d(Im£ n ) . When £„ is centered at zero, then 

«Q) = / ZnPn(tn)dlin = 0. 

The dispersion a n of a distribution p n (£ n ) is given by the equation 

ol EE ((I6J 2 )) = J \UPn{Mn- 

With these definitions, for the stochastic correlation function we have 

((r(W))) = E^<(*)^(*')- 

n 

In the case of white noise, all dispersions <r„ are equal to each other, so that we may write a\ = 2D. Thence, the 
stochastic correlation function becomes 

(<^(0>> = 2£ 

n 

where it is taken into account that is real. From here, because of the completeness of the basis {<p n (t)}, we obtain 
Eq.(59) which was used in calculating expressions (62)-(65). 
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Table Captions 

Table 1. 

The characteristic time t , in seconds, for several parameters £ and different initial conditions. 



Table 2. 

The characteristic time t\ , in seconds, for the atoms initially located at the centre of the trap and for different 
parameters £ . 
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